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We study the spin-transfer torque and magnetoresistance of a ferromagnet| superconductor] ferromagnet spin- 
valve, allowing for an arbitrary magnetization misorientation and treating both s-wave and d-wave symmetries 
of the superconductor. We take fully into account Andreev reflection and also the spin-triplet correlations that 
are generated when the magnetizations are non-collinear. It is found that the torque and magnetoresistance are 
both strongly enhanced when topological zero-energy states are present at the interfaces, which is the case for d- 
wave superconductors with a crystallographic orientation of [1 10] relative to the interface (d^-wave symmetry). 
Moreover, we find that the magnetoresistance displays a strong oscillatory and non-monotonous behavior as a 
function of ds/£ where ds and £ are the inteiiayer width of the superconducting region and the superconducting 
coherence length, respectively. This feature is also attributed to the crossover from layers of size ds ~ 2£ to 
layers of size ds 2£, where the contribution to transport from zero-energy states gradually vanishes. 

PACS numbers: 74.20.Rp, 74.50.+r 



I. INTRODUCTION 

The magnetization dynamics induced by ele ctric curr ents is 
a topic presently benefitting from great interesf^nEHil Typi- 
cally, the dynamics is studied in spin-valve or magnetic tunnel 
junction setups consisting of two or more ferromagnetic lay- 
ers with misaligned magnetization orientations. These layers 
routinely consist of segments with free and fixed magnetiza- 
tion directions separated by non-magnetic spacers in order to 
overcome the exchange coupling between them. While the 
magnetoresistance (MR) of such structures has been the focus 
of most investigations so far, there are other intriguing phe- 
nomena occuring in this type of systems that harbor a great 
potential for applications. One of these phenomena is the 
so-called spin-transfer torque^. The basic concept of spin- 
transfer torque is illustrated by considering a spin-current in- 
cident on a ferromagnetic layer with a magnetization that is 
misaligned compared to the polarization of the spin-current. 
Upon entering the ferromagnetic region, a portion of the in- 
coming spin-current will in general be absorbed by the ferro- 
magnetic order parameter, causing the magnetization to pre- 
cess. This suggests that if the incoming spin-current is large 
enough, it may actually switch the magnetization direction of 
the second ferromagneP, an effect which has been observed 
experimentally^^!. The spin-transfer torque then arises due to 
the non-conservation of the spin-current. 

While the conventional spin-valve structure of type 
ferromagnet | normal metal | ferromagnet (F|N|F) has been the 
subject of much investigation, magnetization dynamics and 
spin-transfer torques in the superconducting (S) analogue 
junction, F|S|F junction, has received less attention so far. 
Only a few works have studied effect s related to the spin- 
transfer torque in such heterostructures^ 12 l 13 l 14 l 15 l In partic- 
ular, the authors of RefP-U calculated the spin-transfer torque 
in an F|S|F spin-valve for an arbitrary misorientation angle 
between the ferromagnets when disregarding the effect of 
Andreev-refiection. However, Andreev-refiection is in general 
present, and should be expected to heavily influence the result- 
ing spin-transfer torque. Also, equal-spin pairing (S z = ±1) 



superconducting proximity amplitudes are induced i n the sys- 
tem whenever the magnetizations are non-collinear! 16 ! 17 !^^ 
Therefore, it remains to be clarified precisely how the super- 
conducting correlations affect the spin-transfer torque in F|S|F 
spin-valves. In addition, to the best of our knowledge there 
exists no study of the spin-transfer torque in F|S|F systems 
where the superconductor has an unconventional pairing sym- 
metry, such as d-wave, which would be relevant for spin- valve 
setups with high-T c cuprate superconductors. 

In this paper, we calculate the spin-transfer torque in a 
F|S|F spin-valve (see Fig. [TJ and fully take into account both 
Andreev reflection and triplet correlations that are present in 
the system. As mentioned, we will allow for an anisotropic 
superconducting order parameter such as d-wave in order to 




FIG. 1: (Color online) A spin-valve setup with misaligned magne- 
tizations in the ferromagnetic regions. We consider both an s-wave 
and d-wave symmetry in the superconducting region. 
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investigate how this influences the spin-transfer torque. More- 
over, we investigate the MR effect of the spin-valve, with par- 
ticular emphasis on the influence of topological zero-energy 
states which may be pre sent when the superconductor has a 
d-wave symmetr y 2 ™ 21 ! 22 !. Our main finding is that both the 
spin-transfer torque and MR are enhanced considerably when 
topological zero-energy states are present at the interfaces, as 
is relevant for the d xy -wave case. Also, the MR exhibits a 
non-monotonous behavior as a function of the width ds of the 
superconductor, an effect which we will show pertains to the 
existence of the zero-energy states. 

We organize this work as follows. In Sec. [TTJ we introduce 
the theoretical framework used to obtain our results, which are 
presented in Sec. [TTTJ We discuss our results in Sec. |IV| and 
give concluding remarks in Sec. [V] 



H. THEORY 

Our model is shown in Fig. [T] We consider transport in 
the x — y plane, while the magnetizations of the F layers are 
allowed to rotate in the y — z plane. This is the most rele- 
vant setup experimentally. We are considering a spin-valve 
setup where an electrical current flows from ferromagnetic re- 
gion Fl (x < 0) where it becomes spin-polarized. The current 
then passes through a superconducting region (0 < x < ds) 
and into a second ferromagnetic region F R (x > dg), on 
which it exerts a spin-torque. The Bogolioubov-de Gennes 
(BdG) equations relevant for our system are formally derived 
by means of diagonalizing the Hamiltonian. In their full form, 
the BdG equations read: 



Hq A(x)ia2 

[A(x)]*lff 2 



HI 



Ip = Elp, 



(1) 



where we have defined the normal-state Hamiltonian Hq and 
wavefunction ip: 



V 2 



2m(x) 



/x(x) 1 — h{x) ■ <t, ip = 



H 



(2) 



The superconducting order parameter is non-zero only for 
< x < dg and the magnetic exchange field is given by 
h(x) — h-L for x < and h(x) = Tir || z for x > ds- Similar 
considerations apply to the chemical potential fi and the effec- 
tive mass m. We denote a 2 x 2 matrix with . ? ., and the super- 
script "T denotes the matrix transpose. The exchange field is 
roughly related to the magnetization M as M ~ —iibNqH 
in the quasiclassical approximation. The elements u„ and v a 
of the wavefunction may be interpreted as electron- and hole- 
like components with spin a, respectively. We have used units 
fi,=c=l and from now on use a real gauge for the superconduct- 
ing order parameter. 

We will now use an extended Blonder- Tinkham-Klapwijk 
(BTK)-formalism 23 to construct the wavefunctions which ac- 
count for the transport properties of the spin-valve, using re- 
sults from RefP^J to account for the misalignment angle be- 
tween the magnetizations. A similar approach was used in 
RefPto investigate the spin-transfer torque in a conventional 
F|N|F spin-valve. For an incident plane-wave from the F L re- 
gion with spin a with respect to the magnetization of Fl at an 
angle of incidence 9 measured from normal to the interface, 
we obtain the following wavefunction in the left ferromag- 
netic region: 



fcos(ip/2)( S] e lk l x + r\e- lk ' x ) + ism(ip/2)( Sl e lk ? x + rj;e- lk ? x )\ 
ism(tp/2)( S] e lk e x + rle~ lk l x ) + cos((p/2)( Sl e lk ? x + rj;e- lk ? x ) 
r\ cos(ip/2)e lk h x — r^isin((^/2)e lfe i a; 
\ -rlism(ip/2)e lk l x + r l h cos((p/2)e lk t x ) 



(3) 



Here, ip is the misalignment between magnetizations of F lay- 
ers. In the superconductor, one may write: 



*s = 



a ± u ± e ±lq " x + d ± v ± j ± e ±1< > hX 
V I b±u±e ±iq - x - c ± v ± -f±e ±lclhX 

±iqhX 



± 



-b±v±j±e ±lq < lX + c±u±e 
a ± v ±1 * ± e ±lq " x + d ± u ± e ±lqhX 



(4) 



whereas the wavefunction in the right ferromagnet reads: 



Fr 



4e lk - x 
t\e- lk l x 



(5) 



Above, {sj = 1,S| = 0} for a =\ and {s| = 0, S| = 1} 
for a =\. In the superconducting regions, the wavefunctions 
are superpositions of left- and right-going electron-like and 
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hole-like quasiparticles, with wave-vectors 

q e [h] = ^2m s ( Ms + [-Ve 2 - |A(0 s )H ^fc2 . (6) 

The momentum perpendicular to the interface is conserved, 
and given as 



\/2mp((ip + ah + e) sin 



(7) 



for incoming particles with spin a. The wavevectors in the 
ferromagnetic region are: 



K[h] = j2m F (^ F + ah+ [-]e) - ft* . 



(8) 



Note that we are distinguishing the effective electron mass and 
the Fermi level in the ferromagnetic and superconducting re- 
gion, in order to employ realistic values for these parameters 
in what follows. The wavevectors g e r/,i in the superconduct- 
ing region are written down under the quite general assump- 
tion that the gap satisfies |A(fc x , k y )\ — |A(— k x) k y )\, which 
is valid for both s-wave, d^-wave, and d x 2_ y 2 -wave sym- 
metries. The superconducting coherence factors are therefore 
given by 



u± [v±] 



1 

T2 



i + H 



Ve 2 -\A(e s )\i 



1/2 



(9) 



where the angle of propagation in the superconducting region 
is given by 



Os 



i(*x/\/2ms/is). 



(10) 



The transport properties of the spin-valve are expected to be 
sensitive to the internal phase of the superconducting order 
parameter, which varies with 8 in the d-wave case. To capture 
this effect, we have introduced the phase-sensitive factors 



7± = A(0f )/|A(0f )|, 0+ = e s ,6s=7r- 9 S . 



(11) 



For the s-wave symmetry we set A(9) = Ao, while for the 
rf-wave case we set A(9) = Ao cos[2(# — a)] with a = for 
d a! 2_ 1/ 2-wave and a = ir/4 for d^-wave pairing. 

The presence of a ferromagnetic material may lead to a 
spin-dependent barrier potential at the interface. In general, 
there may also be spin-flip scattering processes at the inter- 
face due to magnetic inhomogeneities or misaligned moments 
in the interface region 26 . Such spin-flip scattering processes 
may induce equal-spin superconducting spin-triplet correla- 
tions. Since this type of correlations appear anyway due to 
the possibility of a misalignment between the magnetizations 
in the two F layers, we will not include such spin-flip pro- 
cesses at the interface. Quantitatively, the influence of interfa- 
cial spin-flip processes compared to non-collinear magnetiza- 
tions may be different. Qualitatively, the effects brought about 
by the two should however remain the same. Here, we will 
focus primarily on the effect of the misalignment <p between 
the F layers, and we therefore do not write any explicit spin- 
dependence of the barrier potential. Note that the scattered 
particles nevertheless experience spin-dependent phase-shifts 
due to the exchange potential in the F layersPS 



The scattering coefficients in Eq. |3]l are determined by 
using appropriate boundary conditions. Assuming a barrier 
potential of Vq at each of the F|S interfaces, we may define a 
dimensionless measure Z of the transparency: 



Z = 2mFVo/ y / 2mpfj,F- 
The boundary conditions now read: 



(12) 



# R = *s, <9x(# R - *s) = 2m F V ipF R . (13) 

The scattering coefficients are obtained numerically by setting 
up a system of equations 



Ax = b, 



(14) 



where the matrix A and the vector b are constructed by insert- 
ing Eqs. <j3j, (4) and (5) into Eq. ( 13 1. Above, the vector x 
contains all the scattering coefficients. Once these have been 
determined, the wavefunctions are known everywhere in the 
system. With these in hand, we are in a position to calculate 
both the spin-transfer torque and the MR of our system. Let 
us investigate the former in some detail to begin with. The 
spin-operator may be defined as 



S = V t diag( -, < 7*)V/2, 



(15) 



and satisfies the following continuity equation in the super- 
conducting region: 



as 

~dt 



0, 



3s 



2m, 



-Im{*sV[diag(<r,flr*)* s ]}. 



(16) 



Note that there is no source term on the right hand side of the 
continuity equation, in contrast to the corresponding continu- 
ity equation in the ferromagnetic regions. The spin-current 
density is not conserved in the ferromagnetic regions since it 
is partially absorbed by the ferromagnetic order parameter as 
a spin-transfer torque. In the superconducting region, how- 
ever, the spin-current density is conserved and is in general 
a tensor, since it has both a direction of flow in real space 
and a direction in spin space. When the incident spin-current 
from the superconducting region impinges on the interface to 
the right ferromagnet Fr, its transverse component must be 
absorbed since the spin-current deep inside the ferromagnetic 
region must be polarized parallell to its magnetization. There- 
fore, the spin-transfer torque density is equal to the compo- 
nents of the spin-current in the superconductor perpendicular 
to the i-direction, since Mr \\ z. The spin-transfer torque 
may be further decomposed into an in-plane and out-of -plane 
component, rn and t±, with respect to the plane with normal 
vector 



n = (M L x Mr)/ | M L x M r | 



(17) 



spanned by the magnetizations. For our geometry, one in- 
fers that the in-plane component is simply the y-component of 
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the spin-current in the superconductor while the out-of-plane 
component is the A-component. Therefore, we may calculate 
the torque in our system by first obtaining the scattering co- 
efficients, constructing the spin-current of Eq. ( 16 1, and then 
using 



3s -V, r J- =3s 



(18) 



Note that for tunnel junctions with low interface transparen- 
cies, the absorption of the transverse components of the spin- 
current is much less efficient than for metallic interfaces. 
In that case, the spin-current may retain components that 
correspond to a precession around the magnetization axis 
for a larger distance inside the ferromagnet than in high- 
transparency case. Here, we shall assume that the ferromag- 
netic reservoir is large enough to absorbs the transverse com- 
ponents fully. 

In order to gain insight into how the spin-transfer torque is 
affected by superconducting correlations, and in particular the 
existence of topological zero-energy states in the d xy -wave 
case, we proceed with an approach similar to that of RefP. 
We consider incident electrons at the Fermi level (e = 0) 
from the left ferromagnet, adapting a free-electron description 
for these conduction electrons which allows us to employ Eq. 
Q. In order to properly account for classical dephasing, we 
take into account the full range of transverse wavevector com- 
ponents. The spin-transfer torque on the right ferromagnetic 
region is then obtained by computing the spin-current in the 
superconductor for both incident spin-| and spin- j electrons 
at Fermi level, and adding their contribution with a weight- 
factor of 



P CT = (1 + <rh/n F )/2 



(19) 



meant to reflect the different magnitude of the density of states 
due to the finite polarization of the ferromagnet. 

The contribution to the spin-current from a given angle of 
incidence gives rise to an oscillatory behavior in the compo- 
nents of the spin-current which are perpendicular to the mag- 
netization. In our case, these are the x- and y-components 
since Mr \\ z. This oscillatory behavior represents a spin- 
precession around the z-axis upon penetration into the fer- 
romagnetic region with a period of precession which is very 
short for transition metals (a few atomic lattice spacings) 6 . 
However, this oscillatory motion dies out due to classical de- 
phasing which occurs when taking into account all angles 
of incidence of the spin-current. The point is that electrons 
reaching a given distance inside the magnet must have trav- 
elled different path lengths to get there, leading to a destruc- 
tive interference between the phases. 

In addition to considering the spin-transfer torque, we also 
calculate the MR of the superconducting spin-valve depicted 
in Fig. [T] At zero temperature, it may be defined as 



MR: 



Gp — G p. 



(20) 



where Gp/ap is the conductance in the parallell/antiparallell 
(P/AP) magnetization configuration. Whereas the MR has 



been studied for superconducting spin-valves with conven- 
tional s-wave symmetries in several previous works, the case 
of a d-wave superconductor sandwiched between two fer- 
romagnets has received little attention 29 ! 3 ™ 3 ' I J n comparison. 
Here, we will show that the unconventional d-wave symmetry 
can induce a strongly enhanced MR effect, in particular when 
there is a formation of zero-en ergy surface states. Following 
the calculations of Refs ! 32 l 33 l 34 l 35 [ we express the normalized 
conductance of our system when a voltage eV is symmetri- 
cally applied to the junction with eV/2 in Fl and (— eV/2) in 
F R as: 

i.jr/2 



G 



E 

Act 



6.6 cos 6P„ 



-7r/2 



Re{^}, 



,A|2 



MMi|fA|2 



(21) 



where k„ 



is the a; -component of the incident wavevector. It 

that the contri- 



should be implicitly understood in Eq. (21 



bution from incident particles with both spin-| and spin- J. are 
taken into account, leading to different scattering coefficients 
and (the index a has been suppressed in these quantities 
for the sake of notation). Previous theoretical works have con- 
sidered the MR in conventional (s-wave) F|S |F spin-valves for 
either exclusively P or AP configurations 3 - 3 -! 3 -^, while recently 
RefPSl considered an arbitrary magnetization configuration. 
Here, we will investigate the influence of a d-wave supercon- 
ducting order parameter on the MR, with particular emphasis 
on the role of zero-energy states. While it was shown in Refill 
that the presence of zero-energy states in d-wave supercon- 
ducting spin-valves lead to an anomalous voltage dependence 
of the MR effect, it remains to be clarified how the barrier 
transparency and the width of the superconducting layer affect 
the MR. The width of the layer can be controlled accurately 
in the fabrication of experimental samples, while the interface 
transparency may be expected to vary from sample to sam- 
ple. The spin injection properties of heterostructures with un- 
conventional superconductors (p-wave and d-wave) have also 
been studied in Refs! 36 | 37 | 38 [ 



III. RESULTS 

We now discuss our choice of parameters. The ferromag- 
netic regions will be modelled by taking /.ip = 3 eV and 
h = 0.9 eV, which should be appropriate for a transition metal 
ferromagnet. Note that the energy-scale of this exchange field 
(~ 10 4 K) is much larger than room-temperature (~ 300 K). 
In the superconducting region, the Fermi level and gap magni- 
tude depend on whether we have a type I s-wave superconduc- 
tor or a type II d-wave superconductor. For a standard BCS 
superconductor such as Al, we may set A = 1 meV, vp = 10 6 
m/s, and [1$ = 10 eV. On the other hand, for a system such as 
a high-T c cuprate superconductor (e.g. YBCO) featuring a d- 
wave order parameter, a representative choice would A = 20 
meV, i>i? = 5xl0 4 m/s, and fig = 1 eV. We use the for- 
mer set of parameter for the s-wave case, and the latter set of 
parameters for the d xy - and d x 2_ y 2 -cases. Moreover, we will 
consider the ballistic limit, and employ a thickness of the su- 
perconducting region equivalent to several coherence lengths. 
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FIG. 2: (Color online) Plot of the spin-transfer torque and its dependence on the misorientation angle <p. We give results for three different 
symmetries of the superconducting order parameter with interfaces characterized by the barrier parameter Z. Here, we have fixed ds/£ = 3.0. 



For the above parameters, we find that £ ~ 1 nm in the d- wave 
case, which indicates that the thickness of the superconduct- 
ing layer should be of several nm, which is experimentally 
feasible. The electron mass is taken to be close to its bare 
value in the F regions, mp = 0.5 MeV, while in the supercon- 
ductor it is calculated according to the free-electron dispersion 
uig = 2/j,g/r%, where r v — X/vp is the ratio of the velocity 
of light and the Fermi velocity (note that we set c = 1). 

In order to justify a non-selfconsistent approach (see also 
discussion in Sec. |IV| >, we consider for the spin-transfer 
torque a rather thick layer ds/£ = 3.0 and consider several 
values for the barrier transparency Z. Also, we focus only on 
small bias voltages eV/ Ao <C 1 when calculating the MR, 
since larger voltages eV ~ Ao lead to a substantial spin ac- 
cumulation in the superconducting region which strongly sup- 
presses the gap and finally destroys superconductivity itself 39 . 
For eV/Ao <C 1, a self-consistent solution of the gap should 
not be requirecPD. 

A. Spin-transfer torque 

We now proceed to investigate the behavior of the spin- 
transfer torque in the spin-valve setup. The normalization 
constant 

r = v / 2m F fi F /m s (22) 

is introduced below. First, we plot the magnitude t of the 
spin-transfer torque versus the misorientation angle <p in Fig. 
[2] The results are shown for several order parameter symme- 
tries in the superconducting region. The s-wave and d x i_ y i- 
wave symmetries do not feature zero-energy states whereas 
the d xy -wa\e symmetry does. From Fig. [2] one infers two im- 
portant pieces of information. Firstly, it is seen that the depen- 
dence on the misorientation angle essentially follows the usual 
sinusoidal behavior, regardless of the symmetry of the super- 
conducting order parameter. Secondly, whereas the magni- 
tude of the spin-transfer torque rapidly decreases with increas- 
ing barrier strength Z in the s-wave and d x 2_ y 2-wa\e cases, 
the torque is strongly enhanced and more resilient towards an 



increase of Z in the d^-wave case. In fact, in this case the 
spin-transfer torque shows a weak maximum as a function of 
Z before it slowly decreases for large values of Z. This is 
a direct consequence of the existence of zero-energy surface 
states, as will be discussed below. 

To investigate this matter further, we plot in Fig. [3] the 
spin-transfer torque as a function of Z with a fixed misori- 
entation angle <fi = ir/2. Since the magnitude of the torque 
differs greatly for the various symmetries, we use a logarith- 
mic scale for the ordinate axis. As seen, the torque has a 
smooth peak at Z ~ 2 for the d xy -wavs symmetry. Also, 
its magnitude is considerably larger - the ratio of the MR in 
the d xy -wave and s-wave case is ~ 10 3 for Z ~ 2. In con- 
trast, the torque decays monotonically with Z for all values of 
Z > in the s-wave and d x 2_ y 2-wa\e cases. These symme- 
tries do not host zero-energy states at the interfaces, and thus 
an increase in barrier strength simply amounts to a reduction 
of the current since tunneling becomes prohibited, and con- 
comitantly the magnitude of the torque decreases. However, 
when zero-energy states are present the quasiparticle current 
through the superconductor is strongly enhanced at the Fermi 
level compared to e.g. the fully gapped s-wave case, and the 
resulting spin-transfer torque increases in magnitude. This ex- 
plains the monotonic decay of the spin-transfer torque in the 
s- and d x 2_ y 2 -cases when Z is increased, as well as the large 
difference in overall amplitude of the spin-transfer torque be- 
tween, on the one hand the s- and d x 2_ y 2 -cases, and on the 
other hand the d^^-case. The nonmonotonic variation of t/to 
as a function of Z for the d XJ( -case is a more subtle matter, to 
which we will return in Sec. |IV] 

Let us also address the decay of the spin-current inside the 
ferromagnetic region. As previously mentioned, the trans- 
verse components of the spin-current will in general precess 
around the magnetic order parameter, but decay rapidly once 
inside the ferromagnetic region. In Fig. |4] we investigate how 
this decay depends on the barrier parameter Z for various or- 
der parameter symmetries in the superconducting regions. We 
only give the results for the a; -component of the spin-current, 
as the result for the y-component is basically the same. Note 
that the ^-component of the spin-current remains constant in 
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the ferromagnetic region. As seen, the presence of super- 
conducting correlations lead to no qualitatively new behavior 
with respect to the spatial decay of the transverse components 
of the spin current. The spin-current oscillates rapidly on a 
scale comparable to several Fermi wavelengths l/kp where 
we have defined kp = ^/2mp/j,p. In terms of magnitude, it 
is seen that the spin-current is strongly enhanced in the d xy - 
wave case compared to the s-wave and d x 2_ y 2-wave symme- 
try. 



B. Magnetoresistance 

Next, we consider the MR effect in our superconducting 
spin-valve setup, defined by Eq. (20i. In Fig. [5] we plot the 
MR in percentage versus applied bias voltage in the regime 
eV <C Ao where one may disregard spin accumulation. Sim- 
ilarly to the spin-transfer torque, the d xy -wave case behaves 
qualitatively differently from the s-wave and d x 2_ y 2-wave 
symmetries. In the former case, the MR decreases in magni- 
tude upon increasing the bias voltage, while in the latter case 
the MR increases slowly. This may again be understood by 
the enhancement of the quasiparticle current at eV — due to 
the presence of zero-energy states. When eV — > 0, it is seen 
that the MR is enhanced by a factor ~ 10 3 in the d xy -wave 



case compared to the fully gapped s-wave case, similarly to 
the increase of spin-transfer torque discussed in the previous 
section. 

Interestingly, we find that the MR shows a non-monotonic 
behavior upon changing the interlayer width d$ of the su- 
perconducting region in the d xy -wave case. This is shown 
in Fig. [6] where we plot the MR as a function of ds/S, for 
eV/A = 0.01. In the case Z = (perfect barrier trans- 
parency), the MR decays just like in the s-wave case, only 
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FIG. 3: (Color online) Plot of the spin-transfer torque and its depen- 
dence on the interface transparency Z. Here, we have set ds /£ = 3.0 
and if = 7r/2. 



with a much enhanced magnitude. However, for Z ^ the 
MR undergoes a smooth maximum at a finite value of ds /£ , 
before the decaying behavior sets in. 

We propose the following explanation for this observation. 
The topological zero-energy surface states are bound to the 
surface, and penetrate a distance comparable to £ into the su- 
perconducting region. Therefore, they cannot contribute to the 
transport properties of spin-valves with widths <i,g/(2£) ^> 1. 
This is the reason for the decay which sets in at a finite value 
of ds in Fig. [6] signalling precisely the loss of contribution 
to transport from such zero-energy states. It is also seen that 
the MR displays rapid oscillations as a function of ds, per- 
taining to the formation of resonant transmission states in the 
superconducting region. We have checked that upon increas- 
ing the voltage eV, the value of ds/£, at which the MR peaks 
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FIG. 4: (Color online) Plot of the decay of the transverse (x- 
component) part of the spin-current in the right ferromagnetic region 
for several values of the interface transparency Z. Here, we have set 
d s /S, = 3.0 and<^ = tt/2. 
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FIG. 5: (Color online) Plot of the MR for several choices of symmetry for the superconducting order parameter and several barrier transparen- 
cies Z. Here, we have set ds/£ = 3.0. 



is shifted towards smaller values, which is consistent with the 
fact that the influence of the zero-energy states gradually van- 
ishes when moving away from the Fermi level. 

Note that the expression for the conductance in Eq. (21 1 
applies to a situation where a voltage is symmetrically ap- 
plied to the F|S|F junction, with incoming electrons/holes 
from the left/right side, and corresponds to a situation where 
a current is driven through the entire system. This is distinct 
from a situation where a positive bias eV is applied to one 
of the F leads while the superconductor and the other lead 
are both grounded. In the present case, we are interested in 
how the spin-transfer torque and MR are altered depending on 
the magnetization configuration and pairing symmetry in the 
superconductor. This is probed simply by passing a current 
through the F|S|F junction. If the main aim was to investigate 
crossed Andreev reflection signals, it would be more benefi- 
cial to study the aforementioned situation where a voltage is 
applied only between one of the leads and the superconductor 
and then observe a possible current-response in the other lead. 




FIG. 6: (Color online) Plot of the MR for the d xy -wave symmetry 
with eV7A = 0.01. 



Such an effect would be indicative of non-local scattering pro- 
cesses such as crossed Andreev reflection. 



IV. DISCUSSION 

Throughout this study, we have made some simplifying as- 
sumptions in order to proceed with a partially analytical ap- 
proach. For instance, we have adopted the standard step- 
function approximation of the superconducting order param- 
eter Aq at the F|S interfaces, which should be a reasonable 
approximation as long as ds > 2£. A full numerical approach 
would most likely reveal that A is subject to a partial deple- 
tion near the interfaces, but we do not expect that this will 
alter our results qualitatively. Also, we have disregarded any 
thermal effects and their influence on the spin-transfer torque 
since we consider the zero-temperature limit. The presence of 
a superconductor demands that T < T c where T c is the super- 
conducting critical temperature, and it is thus reasonable that 
thermal fluctuations should be neglible in the regime T -C T c . 
In conventional experiments at room- temperature with F|N|F 
spin-valves, thermal fluctuations play an important role in par- 
ticular close to the P or AP alignment ip = {0, 7r} since no 
spin-transfer torque is possible for perfect P or AP alignment. 
However, by considering fluctuations around a magnetization 
axis, a spin-transfer torque becomes applicable. Typically, the 
magnetization in one of the F layers is considered to be static 
with a high magnetization density and correspondingly weak 
fluctuations, while the other one is dynamic and considerably 
influenced by the thermal fluctuations^. 

In Ref. 40 , it was shown that in the idealized situation of 
fully polarized ferromagnets (i.e. half-metals) and perfect in- 
terfaces (Z — 0), the crossed Andreev reflection (CAR) pro- 
cess could lead to very large MR values in an F|S|F spin- valve 
when the superconductor had a conventional s-wave symme- 
try. When the superconductor has a d xy -wave symmetry, we 
find that at normal incidence 8 = the probability for a CAR 
process is exactly zero independent of {Z,ds,eV} and re- 
gardless of whether the magnetization configuration is in the 
P or AP alignment. In the d x 2_ y 2-wa\e case, however, the 
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probability for CAR scattering remains finite at normal inci- 
dence. This finding is in agreement with the results reported 
very recently in RefFH namely that elastic-cotunneling (EC) 
processes are favored along the nodal lines of the supercon- 
ducting order parameter while CAR processes is enhanced 
along directions where the magnitude of the gap is maximal. 
Since = corresponds to the nodal line of a d xy -wa\e order 
parameter, the vanishing CAR amplitude is consistent with the 
above statement. 

The absence of CAR at normal incidence for the d xy -wave 
symmetry suggests that the large MR effect is mainly at- 
tributed to quasiparticle transport through the superconductor 
facilitated by zero-energy states 21 . However, in order to draw 
a firm conclusion about this, we must take into account all an- 
gles of incidence 6. To investigate this quantitatively, let us in- 
troduce some helpful quantities. Consider now only incoming 
spin-| electrons, which dominate the transport properties due 
to the polarization factors P a and the strong exchange field 
h/fXp ~ 0.3. For a given angle of incidence 9, the probability 
of an elastic co-tunneling process (EC) and a CAR process are 
given by 

P Ec{ e) = ^%| 2 , PcarW = ^^N 2 , (23) 

"'inc "'inc 

respectively. As a total measure of the probability for an EC 
and CAR process, we introduce 



o .• . s-wave .„ ,. 

P alignment AP alignment 



tt/2 



Pt = 



j^Y,Pj(e), j — {EC, CAR}, 
9 e=o 



(24) 



where the summation is over all incoming angles and Ng is 
the number of angles summed over (we set Ng — 157 in what 
follows). In this way, Pj £ [0, 1] is a probability measure 
for the EC and CAR processes when taking into account all 
angles of incidence. We plot Pj for eV/Ao = 0.01 and a 
junction of width ds/£ = 3.0 in Fig. |7]in order to understand 
better how they influence the MR effect of the system. We 
have checked that these quantities remain virtually constant 
in the low energy regime eV/Ao -C 1. Consider first the s- 
wave case. As seen, the EC process dominates for small Z, 
regardless of whether the magnetization configuration is P or 
AR Increasing Z (corresponding to tunnel junction), the EC 
and CAR processes become essentially equal in magnitude. 
Since there is so little difference quantitatively between the P 
and AP configuration one expects a small MR effect for the 
s-wave case, which is consistent with our results presented in 
Sec. ImBl 

Consider now the d xy -wavs case. In the P alignment, the 
EC process again dominates the CAR process, but in contrast 
to the s-wave case this now occurs for all Z. The situation 
changes drastically upon reversing the alignment to AP, where 
the contribution from the CAR process becomes dominant at 
large values of Z. The most important point illustrated in 
Fig. [7] for the d xy -wa\e case is that the EC (CAR) process 
is strongly enhanced in the P (AP) alignment compared to the 
AP (P) case. The EC process is facilitated in this context pre- 
cisely due to the existence of zero-energy states which exist 
near the interface and extend a distance of order ~ £ from each 
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FIG. 7: (Color online) Plot of the probability measures for EC and 
CAR processes as a function of the barrier parameter Z, with fixed 
eV/Ao = 0.01 and d s /£ = 3.0. 



interface into the bulk of the superconductor. Simultaneously, 
the CAR process is strongly suppressed in the P alignment be- 
cause of the strong exchange field that makes spin-| holes less 
available in the right ferromagnet. In the AP alignment, how- 
ever, the CAR process can be enhanced by the existence of 
zero-energy states since there are more spin-| holes available, 
with a concomitant reduction of the EC process. The conclu- 
sion is that whereas the EC and CAR processes display a weak 
sensitivity to the magnetization configuration in the s-wave 
case, the presence of zero-energy states in the d^-wave case 
strongly enhances the EC and CAR processes. This enhance- 
ment is not present in the s-wave junction. Upon combining 
the enhancement of EC case with the suppression of CAR in 
the P alignment of the dry-wave, we obtain a very large MR 
effect. Note that it is not the dominance of EC over CAR itself 
that leads to the large MR effect - it is the difference between 
these contributions when comparing the P and AP alignments. 

Finally, note that although the CAR process may dominate 
over the EC process in a certain parameter regime {e.g. for 
large Z in the AP alignment of Fig. [7J, the conductance of the 
junction is always positive. The CAR and EC processes do 
contribute to the conductance with opposite signs, as may be 



seen by expressing the conductance in Eq. (21 1 with the CAR 
and normal reflection amplitudes via the conservation of prob- 
ability currents, but the conductance still remains positive in 
our case of a symmetrically applied voltage to the junction. 
This follows from simple current-conservation upon apply- 
ing a bias voltage difference between the leads. In a scenario 
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where a voltage is applied only to, say, the left lead while 
the superconductor and the right lead are grounded, the cross- 
conductance may take any sign depending on whether CAR 
or EC dominates the non-local transport properties.^ 

The results of Fig. [7]are also relevant to the non-monotonic 
behavior of the spin-torque transfer as a function of Z found 
in Fig. [3] as we discuss below. The EC process contributes 
positively to the magnitude of the current, hence also to the 
torque, and is seen to be maximal at a finite value of Z in the 
dcy-wave case. 

We propose the following explanation for the reason why 
Pec(Z) is non-monotonic in the d xy -case for parallell align- 
ment, while it is monotonically decaying in the other cases: 
When there are no ZES, the EC process cannot be facilitated 
by any subgap states in the superconducting region and it sim- 
ply decays with increasing Z. Consider now the case with 
ZES, which facilitates the EC process. In the AP alignment, 
the EC process is hindered by the density of states mismatch in 
the two ferromagnets as compared to the P alignment, where 
there is no Fermi-vector mismatch associated with tunneling 
of a spin-er electron from the left ferromagnet to the right fer- 
romagnet via the superconductor. In light of this, it is reason- 
able to expect that the ideal circumstances for the EC process 
is the P alignment in the presence of ZES, which is seen to be 
the case from Fig. [7] Note that the gap nodes coincide with 
normal incidence in the d xy -wave case, while the gap ampli- 
tude is maximal at normal incidence in the d x 2_ y 2-wa\e case, 
which also contributes to the suppression of EC in the latter 
case. The reason for why we only see a non-monotonous be- 
havior of P E c in the P alignment for the e^-wave case is thus 
that it is only in this case that the EC process is actually fa- 
vorable without any suppressing mechanisms except for the 
barrier strength Z. In all the other cases, there are additional 
mechanisms which hinder the EC process, which evidently 
leads to a pure decay of the probability for this scattering pro- 
cess. Although the above discussion has been for the P case, 
it bears upon the non-collinear situation <p ^ {0, it}. The rea- 
son for this is that such a scenario may be considered as a 
superposition of the P and AP alignments, such that the non- 



monotonous behavior should be seen for any ip which is not 
too close to 7r. The exact microscopic mechanisms behind the 
non-monotonous behavior of the EC process are very diffi- 
cult to extract since our approach is partly numerical, but we 
believe the general arguments given above should provide a 
basic understanding of the situation. 



V. SUMMARY 

To summarize, we have investigated the spin- 
transfer torque and the magnetoresistance of a ballistic 
ferromagnet I superconductor I ferromagnet spin-valve, al- 
lowing for an arbitrary magnetization misorientation and 
also including the possibility of a d-wave symmetry for the 
superconducting order parameter. Our approach accounts 
properly for both Andreev reflection and proximity-induced 
triplet correlations that are generated when the magnetiza- 
tions are non-collinear. Our main finding is that the torque 
and magnetoresistance are both strongly enhanced when 
topological zero-energy states are present at the interfaces, 
which is relevant for the d xy -wa\e case. Also, we have 
found a strong oscillatory and non-monotonous behavior of 
the magnetoresistance as a function of the width d$ of the 
superconductor. The reason for this is a gradual vanishing of 
the contribution to transport from zero-energy states as the 
layer size increases from ds ~ 2£ to ds 3> 2£, where £ is the 
superconducting coherence length. 
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